DFT study of nitrided zeolites: Mechanism of nitrogen substitution in HY and silicalite

0021-9517/$ see front matter 2009 Elsevier Inc. A doi:10.1016/j.jcat.2009.10.015 * Corresponding author. Address: Goesmann 222, Amherst, MA 01003-9336, USA. Fax: +1 413 545 4490 E-mail addresses: [email protected] (V. Ag (G.W. Huber), [email protected] (W.C. Conner edu (S.M. Auerbach). We have performed embedded-cluster calculations using density functional theory to investigate mechanisms of nitrogen substitution (nitridation) in HY and silicalite zeolites. We consider nitridation as replacing Si–O–Si and Si–OH–Al linkages with Si–NH–Si and Si–NH2–Al, respectively. We predict that nitridation is much less endothermic in HY (29 kJ/mol) than in silicalite (132 kJ/mol), indicating the possibility of higher nitridation yields in HY. To reveal mechanistic details, we have combined for the first time the nudged elastic band method of finding elusive transition states, with the ONIOM method of treating embedded quantum clusters. We predict that nitridation of silicalite proceeds via a planar intermediate involving a ring with pentavalent Si, whereas nitridation of HY is found to proceed via an intermediate similar to physisorbed ammonia. B3LYP/6-311G(d,p) calculations give an overall barrier for silicalite nitridation of 343 kJ/mol, while that in HY is 359 kJ/mol. Although the overall nitridation barriers are relatively high, requiring high temperatures for substitution, the overall barriers for the reverse processes are also high. As such, we predict that once these catalysts are made, they remain relatively stable. 2009 Elsevier Inc. All rights reserved.